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We report contributions to the nucleon spin and momentum from light quarks calculated using 
dynamical domain wall fermions with pion masses down to 300 MeV and fine lattice spacing 
a ~ 0.084 fm. Albeit without disconnected diagrams, we observe that spin and orbital angular 
momenta of both u and d quarks are opposite, almost canceling in the case of the d quark, which 
agrees with previous calculations using a mixed quark action. We also present the full momentum 
dependence of « = 2 generalized form factors showing little variation with the pion mass. 
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1. Introduction 



Nucleon structure calculations on a lattice are important for both testing QCD as the fun- 
damental theory of quarks and gluons and making predictions and complementing experimental 
efforts that aim to measure the full three-dimensional picture of the proton and nucleon, and to 
explore the origin of the nucleon spin. Lattice calculations of nucleon structure are pursued by 
many groups, see [1, 2, 3, 4] for recent advances. In this work, we report results for quark momen- 
tum fraction, spin and orbital angular momentum from calculations with Nf = 2-\-l domain wall 
fermions. We also report previous Nf = 2+1 calculations using a mixed (domain wall valence and 
Asqtad sea) quark action [5] with substantially increased statistics and reduced enmbars [10]. 



Table 1: Lattice ensembles used in this work. Uncertainties in the pion masses are dominated by the scale 



uncertainty. 



X Lt 


m,i [MeV] 




Sep. T/a 




amji 


# confs 


# meas 




Mixed action {a = 0.1240(25) fm, am^tmnge = 0.050) 






20^ X 64 


293(6) 


3.68 


9 


0.007 0.0081 


0.1842(7) 


463 


3704 


20^ X 64 


356(7) 


4.48 


9 


0.010 0.0138 


0.2238(5) 


631 


5048 


28^ X 64 


356(7) 


6.27 


9 


0.010 0.0138 


0.2238(5) 


274 


2192 


20^ X 64 


495(10) 


6.23 


9 


0.020 0.0313 


0.3113(4) 


486 


3888 


20^ X 64 


597(12) 


7.50 


9 


0.030 0.0478 


0.3752(5) 


563 


4504 




Unitary domain wall {a = 0.0840(14) fm, am^trange = 0.030) 




32^ X 64 


297(5) 


4.06 


12 


0.004 0.004 


0.1268(3) 


616 


4928 


32^ X 64 


355(6) 


4.86 


12 


0.006 0.006 


0.1519(3) 


883 


7064 


32^ X 64 


403(7) 


5.52 


12 


0.008 0.008 


0.1724(3) 


528 


4224 



We perform calculations that span different values of lattice spacing, volume and pion masses, 
see Tab. 1. With the domain wall fermion action, we perform high-statistics calculations using 
300 MeV < < 400 MeV, and with the mixed action using 300 MeV < < 600 MeV. 

Quark momentum and angular momentum can be extracted from the matrix elements of the 
quark energy-momentum tensor between nucleon states: 

{P'\q[Yi^ iD^^ - (trace)]^!/') = U{P') \A\,{Q^)y^^ P^"^ +Bl,{Q')-i-pi^o^}-qa 

^ (1.1) 

+ Cl,{Q')--q^^q^}-{tmce)]u{P), 
Mn J 

where q = P' — P, P = ^{P + P'), = —q^ and A20, B20 and C2 are momentum transfer-dependent 
n = 2 generalised form factors (GFF). The forward value of A2o(0) gives the quark momentum 
fraction inside the nucleon, while the combination of A20 and B20 gives the quark angular momen- 
tum [6]: 

(^)=A2o(0), /, = ^[A2o(0)+S2o(0)]. (1.2) 

Note that the B20 GFF requires calculations with non-zero momentum transfer. In this work we 
concentrate onn = 2 GFFs, with n = 3 GFFs reported elsewhere [8]. 
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2. Renormalization 

The energy-momentum operator in Eq. ( 1 . 1 ) is scale-dependent and its matrix elements have to 
be renormalized for comparison with experiment. For the domain wall ensembles, we calculate the 
renormalization constants using the Rome-Southampton method [7] following the RI'-MOM pre- 
scription. In this scheme, the in- and out-quarks of the pertinent vertex function have the same off- 
shell 4-momentum p = p' which defines the scale of the renormalization point jj.^ = p^. This scale 
must be in the window Agf^p <^ jJ.^ <^ Af^^ = a^^, i.e. in the perturbative region but small enough 
to avoid significant discretization errors. There are indications that such a window does not exist 
for the hybrid action ensembles [8], and we use the renormalization constants from lattice pertur- 
bative calculations [9]. Since the gauge field in hybrid action calculations was HYP-smeared [10], 
we believe that the 1-loop perturbative calculations are a reasonably good approximation, and in 
Sec. 3, we will see close agreement between calculations using perturbative and non-perturbative 
renormalization. 

1.79 
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1.75 
1.74 
1.73 
1.72 
1.71 
1.7 
1.69 



Figure 1: RI/MOM renormalization constant divided by 3-loop perturbative running (scale-independent). 
The extracted values must be multiplied with factors Zj^jZ^/ and {Z^^ /Z^l )2 cev to get the final renormal- 
ization constants. 

We calculate vertex functions of the quark energy-momentum operator for two hypercubic 
representations t'^'' and T^^' [1 1] and the (local) axial-vector current operator = qxYiiYsq.x- We 
use the latter to exclude the quark field renormalization, employing the fact that the conserved 
axial-vector current s/^,^ of the domain wall fermions, with good precision, should have trivial 
renormalization 1, and calculating the ratio Z^/Za from pion matrix elements separately. 

In Fig. 1 we show the "scale-independent" renormalization constants, i.e., the lattice renormal- 
ization constants divided by the running constant Z'^''(/i) calculated perturbatively up to 3 loops in 
continuum theory [12]. In the region 4 < jLt^ < 14 GeV^, the "scale-independent" ratios vary within 
±1%, demonstrating perfect matching between the lattice and the continuum renormalization fac- 
tors. Assuming that the discretization errors contribute as Z{a) = Z(0) +Z'a^iJL^, we extrapolate to 
/X = and rescale our results to the final MS (2 GeV) scheme. The final renormalization factors for 
the two representations are close to each other: z'*!^^^ = 1 .708(4) and z'^g^^^ '^^^^ = 1 .749(2). 
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3. Results 

To compute nucleon two- and three-point functions we use Gaussian-smeared sources tuned 
in Ref. [13]. Tiie source-sink separation T for both calculations is roughly 1.0. .. 1.1 fm which we 
expect to be sufficient to suppress excited state effects at the precision level we aim at. We extract 
the matrix elements of the energy-momentum tensor using the standai^d ratio/plateau technique [14] 
and the overdetermined analysis [13, 14]. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 2: (Color online) Isovector (left) and isoscalar (right) « = 2 Generalized form factors (GFF) of the 
proton from domain wall calculations. 

In Fig. 2 we show the dependence of GFFs on the momentum transfer. The bands con^espond 
to the dipole fits of each form factor and their statistical variation. Using these fits, we extract the 
forward {Q^ = 0) values of the isovector and isoscalar GFFs, as well as their slopes. 

It is remarkable that there are specific features of the GFFs independent of m^i, e.g. 

iq+'^l > iq-^'l ^ 0, iB^o'^l > \B"+''\ ^ 0, (3.1) 

which are also consistent with previous calculations [5] and laige-Nc counting rules [15]. 

We perform chiral extrapolation of the generalized form factors at 2^ = using baryon chiral 
perturbation theory [16]. While the formulas in Ref. [16] predict both the mj^ and dependence, it 
is not clear if they are applicable for > 0.2 GeV^. We therefore restrict our analysis only to the 
forward values and their derivatives obtained using the dipole extrapolation for — )• 0. In Fig. 3 
we show fitting and extrapolation results for the isovector GFFs. The analysis of the isoscalar GFFs 
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Figure 3: (Color online) Chiral extrapolations of isovector GFFs A20, ^20. C20 (left) and their derivatives 
d/dQ^ (right) at = 0. The ChPT predictions are taken from Ref. [16]. 
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Figure 4: (Color online) Isovector (left) and isoscalar (right) quark momentum fraction in the proton com- 
pared to calculations in Ref. [1,2]. 



is similar. The OFF derivatives help to constrain the fits considerably, although less than the GFFs 
themselves. 

In Fig. 4 we show our results for the quark momentum fraction (x) for both domain wall and 
hybrid action calculations. It is remarkable that the results agree very well, although lattice spacings 
and renormalization procedures are different. In addition we show the comparison of our results 
to those of Ref. [1] where Nf = 2 Clover fermions were used and Ref. [2], which used Nf = 2 + 1 
domain wall fermions with (3 = 0.114 fm. It is interesting that Nf = 2 + 1 domain wall results tend 
to agree at smaller pion masses, while Nf = 2 Clover calculations disagree by « 15%. Our chi- 
rally extrapolated value of (x)"^^ = 0.204(6) significantly overshoots the phenomenological value 
0.155(5). This disagreement can be the result of either unreliability of the chiral extrapolation or 
presence of excited states in the matrix element calculation. At the same time, the isoscalar momen- 
tum fraction (x) agrees well with the phenomenology, albeit without disconnected contractions. 
This fact may be an indication of the smallness of the disconnected contributions. 

In Fig. 5 we show quark angular momentum (disconnected contractions are not included), as 
well as separate contributions of u and d quark spin and orbital angular momentum (0AM) to 
the proton spin. Although not shown, results from the hybrid action calculations [10] agree very 
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well with the data in Fig. 5(right). The calculations with domain wall fermions confirm previous 
qualitative results^ that 



\f'\<.\J"\, 



r u+d I 



<|L"|,|L'^ 



(3.2) 
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Figure 5: (Color online) Left: u and d quark contributions to the nucleon spin from the domain wall and 
hybrid action calculations. Right: u and d quark spin and orbital momentum from the domain wall calcula- 
tions. 



4. Discussion 

We have computed the generalized form factors of the nucleon in fully dynamical QCD with 
Nf = 2 + 1 flavors of chirally symmetric quarks. The pion masses 300—400 MeV are still too heavy 
to be realistic, however, they are arguably close enough for the applicability of chiral perturbation 
theory. Nevertheless, we see disagreement between the extrapolated value of the isovector quark 
momentum fraction (x)"^^ and phenomenology. The presented agreement with the previous cal- 
culations employing different quark action, lattice spacing and operator renormalization procedure 
indicates that the problem may reside in either insufficiency of the order of the chiral perturbation 
theory expansion or presence of nucleon excited state contributions in the correlation functions. 
There are indications for excited state contributions in other calculations as well [17]. 

The disagreement with the Nf = 2 calculations may also indicate that the number of flavors 
plays a crucial role for this quantity, leading to the conclusion that the low energy dynamics of 
QCD may affect (x)"^'' more than other quantities. 
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